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Abstract 

It is known that the Meissner-like effect is seen in a magnetosphere without an electric current in 
black hole spacetime: no non-monopole component of magnetic flux penetrates the event horizon if 
the black hole is extreme. In this paper, in order to see how an electric current affects the Meissner- 
like effect, we study a force- free electromagnetic system in a static and spherically symmetric 
extreme black hole spacetime. By assuming that the rotational angular velocity of the magnetic 
\ field is very small, we construct a perturbative solution for the Grad-Shafranov equation, which 

is the basic equation to determine a stationary, axisymmetric electromagnetic field with a force- 
free electric current. Our perturbation analysis reveals that, if an electric current exists, higher 
multipole components may be superposed upon the monopole component on the event horizon, 
even if the black hole is extreme. 
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I. INTRODUCTION 



It is widely believed that there are supermassive black holes at the centers of galaxies, 
and these are hypothesized to be the central engines for active galactic nuclei (AGNs) and 
gamma ray bursts (GRBs). Two main possibilities are considered as the energy source. One 
is the gravitational energy of accreting matter and the other is the rotational energy of the 
black hole or the accretion disk surrounding it. However, the details of the energy extraction 
process are not clear. It is also not understood well how the energy is converted into that 
of AGNs or GRBs. 

Blandford and Znajek showed that the rotational energy of a rotating black hole can 
be extracted in the form of Poynting flux along magnetic field lines penetrating the event 
horizon [lj], which is known as the Blandford-Znajek (BZ) mechanism. Its efficiency depends 
on the rotational velocity of the black hole and the configuration of the magnetic field: the 
extraction of the rotational energy becomes more efficient the more magnetic field lines pen- 
etrates the event horizon and the more rapidly the black hole rotates. In the BZ mechanism, 
poloidal magnetic fields which penetrate the event horizon play a crucial role for the energy 
extraction as well as for the formation of jets associated with AGNs. In fact, some numerical 
studies reported that Poynting-dominated jets were produced 2|-|4|]. 

Bicak and Janis showed that a magnetic field without an electric current is expelled from 
the event horizon of a maximally rotating black hole j^j. This is analogous to the Meissner 
effect in a superconductor. This effect for a rapidly rotating black hole would decrease 
the efficiency of the BZ mechanism, though the larger rotational velocity of the black hole 
would increase the efficiency. In realistic astrophysical cases, however, there would be plasma 
around the black hole. How the Meissner-like effect is affected by the existence of plasma is 
the main subject of this paper. We clarify the effect of an electric current on the Meissner- 
like effect of an extreme black hole. Komissarov and McKinney studied numerically the 
Meissner-like effect of a Kerr black hole js]. They carried out numerical simulations for a 
highly conductive magnetosphere until it almost reaches steady state, and there was no sign 
of the Meissner-like effect in their numerical results. In this paper, we study how an electric 
current affects the Meissner-like effect by solving a stationary problem analytically. 

Since realistic situations are, in general, very complicated, it is difficult to model them. 
In order to reveal the essence of the plasma effect, we consider a very simple toy model: 
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(i) we consider a stationary, axisymmetric force-free system of the electromagnetic field and 
plasma; (ii) we consider a static spherically symmetric black hole spacetime with a degenerate 
horizon as a background spacetime rather than a rotating black hole. The degenerate horizon 
is the origin of the Meissner-like effect in a vacuum black hole spacetime |7|, and hence, by 
studying the electromagnetic field in this spacetime, we can see whether the Meissner-like 
effect remains even in the case with an electric current. The spacetime considered in this 
paper is known as the Reissner-Nordstrom (RN) spacetime. By these assumptions, the basic 
equations reduce to only one quasi-linear elliptic equation for the magnetic flux function 
called the Grad-Shafranov (GS) equation 8]. 

For the black hole spacetime, the GS equation has three regular singular points: one is at 
the event horizon, and the other two are at the inner and outer light surfaces on which the 
velocities of the magnetic field lines agree with the speed of light. For non-extreme cases, 
one boundary condition is imposed at each regular singular point so that the magnetic field 
is smooth everywhere. However, for a given electric current function, the obtained solution 

_ h 

for the magnetic flux need not be C 1 but at most C 1 [9(. Although numerical C 1 solutions 



have been obtained by iteratively changing the functional form of the electric current 

a mathematically rigorous proof for the existence of a C 1 solution has not yet been presented. 

Furthermore, in the extreme case, two kinds of boundary condition must be imposed at once 

on the event horizon. We shall mention all these difficulties in solving the GS equation in 

§IV. 

As will be shown in §V, the monopole component is a unique configuration of the magnetic 
field on the event horizon if there is not an electric current. Since there is no magnetic 
monopole in nature, this result implies the Meissner-like effect of the extreme RN black 
hole. In order to study the electromagnetic field coupled to an electric current around an 
RN black hole, we use a perturbative method which includes two expansion parameters. One 
of these parameters corresponds to the rotational angular velocity of the magnetic fields. 
Namely, we consider slow-rotating magnetic fields as was first considered by Blandford and 
Znajek The other parameter is the ratio of the distance from the event horizon to the 
horizon radius, since we consider only the vicinity of the event horizon, which includes the 
inner light surface. Although we cannot take into account the outer light surface in our 
perturbative method, we can obtain approximate solutions sufficient to study the Meissner- 
like effect with an electric current. 
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This paper is organized as follows. In §11, we introduce the RN black hole as a back- 
ground geometry. Then we show the GS equation for the RN spacetime in §111; the detailed 
derivation of the GS equation is given in Appendices A and B. The regularity conditions 
for the GS equation and difficulties in solving this equation are described in detail in §IV. 
Using perturbative analyses, we study the cases with and without an electric current in §V 
and VI, respectively. §VII is devoted to summary and discussion. In Appendix C, we show 
the relation between the Kerr-Schild coordinate system and the standard static coordinate 
system of the RN spacetime. In Appendix D, we give a proof of a theorem on the magnetic 
field obtained by the present perturbative method. 

In this paper, we adopt the geometrized units, in which the Newton's gravitational con- 
stant and the speed of light are unity, and the abstract index notation: small Latin indices, 
excluding t and r, indicate the type of tensor, whereas small Greek indices, excluding 9 and 
(f, represent components with respect to the coordinate basis. The exceptional indices t, r, 
8, and ip denote the components of time, and the radial and azimuthal coordinates in the 
spherical polar coordinate system. The signature of the metric is diag[— , +,+,+]. 



II. BACKGROUND GEOMETRY 



We consider a static and spherically symmetric spacetime of the following metric: 

2 

ds 2 = g^dx^dx" = -a 2 dt 2 + T -dr 2 + r 2 (d9 2 + sin 2 9d<p 2 ) (2.1) 

with 

A = (r — r + )(r — r_) and a = , (2-2) 

r 

where we assume r + > r_ > 0. This spacetime is known as the RN spacetime. There are 
two horizons, which are determined by A = 0: r + and r_ represent the radius of the event 
and Cauchy horizons, respectively. The case of r + = r_ := rn is called the extreme case. 



III. GRAD-SHAFRANOV EQUATION 

Maxwell's equations are given by 

V [a F bc] = 0, (3.1) 
V b F ab = 4vrJ a , (3.2) 
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where F ab is the field strength tensor of the electromagnetic field, J a is a current density, and 
Vfe is the covariant derivative 1 . If the field strength tensor is expressed by using a 4- vector 
potential A a as 

F ab = V a A b - V b A a = d a A b - d b A a , (3.3) 

then Eq. (13.1 j) is trivially satisfied, where d a is the ordinary derivative. 

As mentioned in §1, hereafter, we consider the axisymmetric and stationary electromag- 
netic field in the RN spacetime. In order to make the problem simple, we assume that the 
system field satisfies the force-free condition 

F ab J b = 0. (3.4) 

Hereafter, we focus on the system of only Eqs. (13.21) and (13. 4p . The formulation of a force-free 
electrodynamics field in the black hole spacetime was given by Macdonald and Thorne 
Our formulation is based on their work. 

In the case of a stationary, axisymmetric electromagnetic field, we can define the "angular 
velocity" of the magnetic field as 

n F tg F tr 

Q F := — = — (3.5) 

r Otfl £ Tip 

The reason why can be regarded as the angular velocity of the magnetic field is described 
in Appendix |Al From Eqs. (13 .2\\ and (13. 4ft . the GS equation is obtained as 

*» + -L L ,* + £ + Z. = , (3.6) 



where a prime ' represents a derivative with respect to r 

} d_ ( 1 
d6 \^0~dd 



d ( 1 

V£ = sin0^ \^-~ , (3.7) 



U = (D'+^^9^'^', (3.8) 

W = (d e D + - sin 2 O^P-de* d e ^> + 8n 2 r 2 ^-, (3.9) 
V 2 dty J dty 



1 We only use the RN spacetime as a background geometry rather than the Kerr spacetime, that is, a test 
U(l) field on the RN spacetime is considered in our toy model. Bicak and Dvorak studied perturbations 



of the coupled Einstein-Maxwell system on the extreme RN spacetime without an electric current 15 1. 
The Mcissncr-like effect appeared in this case. 
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where I = and $ are the electric current and the magnetic flux through an axisym- 
metric polar cap, which are defined by Eqs. ( IB4j) and ( 1B5j) . respectively, and D is defined 
by 

D := a 2 — r 2 Qp sin 2 6. (3.10) 

The derivation of the GS equation (13. 6ft is given in Appendices lAl and iBl 

In general, the GS equation for a black hole magnetosphere has three regular singular 
points: one is at the event horizon A = 0, and the other two are at the light surfaces 
defined by D = 0. A light surface is a timelike hypersurface on which the rotational speed 
of magnetic field lines are equal to the speed of light. The inner light surface given by a 
function r = r^-iO) has a spacelike section with spherical topology, whereas the outer one, 
r = tls+(6 i ), has that of cylindrical topology. 

In the case of the Kerr spacetime, the Kerr-Schild coordinate system is often adopted, for 
example, in numerical simulations (e.g., j^l, [l6[ since there is no coordinate singularity 
on the event horizon. Therefore, we give the Kerr-Schild coordinate system for the RN 
spacetime in Appendix In Appendix 0> we show that the singular point on the event 
horizon in the GS equation appears even if we adopt Kerr-Schild coordinates. This is also 
true in the case of the Kerr spacetime, though it is not shown in this paper (we will show 
it elsewhere). As long as a stationary magnetic field is considered, the singular point of the 
basic equation will appear on the event horizon, since there is no timelike Killing vector 
field on or inside the event horizon, or in other words, the stationary configuration cannot 
be realized inside the black hole. 

IV. REGULARITY CONDITIONS AND BOUNDARY CONDITIONS 

In this section, we consider only the case of f2p 7^ 0. Then, by virtue of the symmetry of 
the background spacetime, without loss of generality, we may assume Qf > 0. The case of 
flp = will be treated as specific cases later. 

A. Symmetry Axis 9 = 

On the symmetry axis 8 = 0, both \1/ and / should vanish by their definitions. / is a 
function of and thus, /|<ir = o should vanish. 
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B. Event Horizon 



In Appendix [Cj we give the components of F a b in the Kerr-Schild coordinate system 
(T, $, R, 6). From Eq. flUT3|) . we have 

i^e = (4.1) 

where 

1 2r 2 

.M = — 2Mr-r')n F 9 e $ + — /(f). (4.2) 
27r v 7 sin 6/ 

Since the Kerr-Schild coordinate system is non-singular on the event horizon, Fr® must be 
finite there. Thus, we have Ai\ r=r+ = 0. This leads to 

I+- — : .(2Mr -ri)n F sm9d e y = at r = r + . (4.3) 

The above condition corresponds to the horizon boundary condition derived by Znajek for 

It is seen from the GS equation ( 13. 6 p that, in order that \P' and 



the Kerr black hole 



f are finite on the event horizon r = r + , the following condition should be satisfied 

W 

L e V + — = 0. (4.4) 

However, this condition is satisfied if \& satisfies the regularity condition ( 14. 3ft . and thus no 
additional constraint is imposed by this equation. 

In the extreme case, since the equation A = has a double root r± — rn, -M'| r=rH should 
vanish as well as .M| r=rH = so that Frq is finite on the event horizon. These conditions 



I+^-wxedgV = 0, (4.5) 

47T 



dl T , sin^ „ T , dOp 



imply 



_*' + _^ + _ f ^ = 0, (46) 

on the event horizon r = rn- We can see from the GS equation (I3.6P that, in order that 
f , and \&" are finite on the event horizon for the extreme case, not only Eq. ( 14. 4 ft but also 
the following condition should be satisfied, 

L e y + —)=0 at r = r H . (4.7) 

This condition is satisfied if \& satisfies the regularity conditions (14 . 5 [) and (14.6 j) . and thus 
no additional constraint is imposed by this equation. 



z 
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FIG. 1: Schematic diagram of the magnetosphere of a non-extreme RN black hole r + > r_, where 
p = r sin and z = r cos 0. 

C. Light Surface 

As mentioned, the light surfaces are singular points of the GS equation (13. 6p . In the 
extreme case, the radial coordinates of the light surfaces, which are the roots of the equation 
D = in the domain r > th, are given by 

1 ± VI - 4r H fi F sin0 



r = r LS± 



2f2p sin f9 

In order that W and L^\l/ are finite on the light surfaces, U + W/A must vanish there. This 
requirement leads to the following regularity conditions on the light surfaces: 

,dl 2 



2„2' 



V a d a V = -87r 2 r 



at r = r LS ± 



where 



T/a . = r 2 g ab I + '_ ^2 qq^ 



(4.9) 



(4.10) 
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D. Boundary Conditions in the Case of Non-Degenerate Horizons 



For simplicity, in this subsection, we assume that f2p is constant. By virtue of this 
assumption, the conditions (14. 9p become Neumann boundary conditions on the light surfaces. 

Here, we consider the non-extreme case r + > r_. Let us assume that the functional 
form of has already been determined before solving the GS equation. Then, by solving 
the horizon regularity condition (14. 3p . we obtain ^ on the event horizon. Since ^ is the 
magnetic flux through the polar cap (see Appendix B), \l/ should vanish at 9 = 0. Thus, 
there seems to be no freedom for setting a boundary condition for \1/ in solving Eq. (14.3)) . but 
this is not true. Since 9 = is a regular singular point of Eq. (I4.3|) . there still remains one 
degree of freedom for choosing a boundary value of the second-order derivative of Hence, 
a Dirichlet boundary condition for the GS equation is determined on the event horizon by 
the regularity condition (14. 3p . By imposing this Dirichlet boundary condition at r = r + 
and further Neumann boundary conditions at r = tls- (14.91) and on the equatorial plane 
(e.g., the reflection-symmetric boundary condition dg^/\g =n / 2 = 0), a solution for the GS 
equation (I3.6p is uniquely determined in the domain r + < r < r^s_. By imposing two 
Neumann boundary conditions (14. 9 p on the two light surfaces r = r^s± and the equatorial 
plane 9 = n/2, and a further Dirichlet boundary condition f = on the symmetry axis 
9 — 0, a solution for the GS equation is uniquely determined in the domain r^s- < r < tls+- 
For the domain Tls+ < r < oo, if we impose a boundary condition at \1/ for r — > oo, then we 
can obtain a solution for the GS equation. 

The GS equation can be solved for these three domains, r + < r < r LS _, r L s- < r < r LS+ , 
and r LS+ < r < oo, independently by the above procedure (see Fig.[T|). Thus, for an arbitrary 
electric current the obtained solution for \1/ is, in general, not C l but at most C l ~ 

at the boundaries r = tls±- It was first reported by Contopoulos, Kazanas, and Fendt 
(CKF) that the continuity of the first-order derivative of \1/ at the light surface as well as the 
continuity of \1/ itself determines the functional form of the electric current in the case of 
the pulsar magnetosphere js]. They numerically obtained C 1 solutions for \1/ by an iterative 
method in which both \1/ and the functional form of /(^) are determined simultaneously. 
The CKF method was used by several authors to study the pulsar magnetosphere and they 



showed that this method also suitable for their case of interest 
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141 ]. However, it 



should be noted that a mathematically rigorous proof for the existence of the C 1 solution 
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for the GS equation has not yet been given. 

Although there is at most one light surface in the case of the pulsar magnetosphere, there 
can be two light surfaces if Qp is a non-vanishing constant in the case of the black hole 
magnetosphere (see Fig. [[]). Thus, if we chose the functional form of such that \& is C 1 
in the domain r + < r < Tls+, imposing an asymptotic boundary condition for r — > oo does 
not guarantee the continuity of the derivative of ^ at the outer light surface r = tls+- Thus, 
in order to obtain a solution which is C 1 in the domain r + < r < oo, we need to solve the 
GS equation for the outermost domain tls+ < t < oo as a Cauchy problem with a boundary 
data for ^ and the derivatives of \& on r = rLg + . This implies that we cannot impose the 
asymptotic boundary condition for r — > oo. In general, it is difficult to solve an elliptic-type 
differential equation numerically, such as the GS equation, as a Cauchy problem due to the 
numerical instability. Thus, in the case with two light surfaces, it is difficult to numerically 
obtain a solution for the GS equation in the outermost domain tls+ < r < oo. However, this 
might not be a serious problem, since we may understand whether the Blandford-Znajek 
mechanism works by studying only the domain r + < r < 7"ls+- 

Uzdensky applied the CKF method to the magnetospheres of the Schwarzschild black 



hole [llj and of the Kerr black hole [12|, but he focused on only the cases in which there 
is only one light surface by virtue of a particular assumption on Qp: Uzdensky assumed 
that Qp asymptotically decreases, and hence there is no outer light surface. Thus, Uzdensky 
succeeded in numerically obtaining global solutions without solving Cauchy problems for 
the GS equation. 

E. Boundary Conditions in the Case of a Degenerate Horizon 

Here, we consider the extreme case r + = r_ = rn, which is the main case of interest in 
this paper. In this subsection, we also assume that flp is constant. The horizon regularity 
conditions (14. 5[) and (14.61) give boundary values of ^ and the derivative of Thus, in the 



extreme case, we must solve the GS equation (13. 6p as a Cauchy problem even for the domain 
r H < r < r L s- (see Fig. [2]). As mentioned, it is difficult to numerically solve the GS equation 
as a Cauchy problem, and hence, it seems to be difficult to numerically obtain a solution 
for the physically important domain < r < Tls+ in the extreme case. Further, even 
if we find a procedure for numerically solving the GS equation as a Cauchy problem, the 
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Inner Light Surface r - r L s- 
No boundary condition can be imposed. g 



f Outer 
' Light Surface 
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FIG. 2: The same as Fig. [TJ but for an extreme RN black hole r + = r_ = rn- 



regularity condition on the inner light surface r = tls- may not be satisfied for an arbitrary 
functional form of electric current we must assume the functional form of to 

solve the GS equation as a Cauchy problem, but, in general, the assumed electric current 
does not satisfy the regularity condition on the inner light surface. As a result, it seems 
to be impossible to obtain a solution for the GS equation numerically, which is finite on the 
inner light surface, in the extreme case. In this sense, the perturbative analytic approach 
discussed in §VI is very important. 

We should note that even if we find analytically the electric current / which guarantees 
the finiteness of \l/ and its derivative on the inner light surface, such a electric current / might 
not guarantee the finiteness of both \1/ and its derivative on the outer light surface. This 
implies that either the force-free condition should break down near the outer light surface 
or the rotational velocity should decay far from the black hole so that the outer light surface 
does not exist, as in the situation studied by Uzdensky. 

V. VACUUM CASE 

In this section, we consider the vacuum case 1 = 0. 
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A. rip = on the horizon 



Here, we consider the case of f2p = on the event horizon, r = r + . Even if Op does not 
vanish except on the event horizon, it satisfies 

d e n F = ^-d e ty = at r = r+. (5.1) 

Thus, on the event horizon, dgty = or dQp/d^/ = should hold. In the former case, ^ 
vanishes on the event horizon, implying that the magnetic flux does not penetrate the event 
horizon. In the latter case, from Eq. (I3.6p . in order that and ty" are finite on the 

event horizon, the following equation should be satisfied: 



AY 



L e * + rV^— J =0 at r = r + . (5.2) 

In the extreme case, since A' = also holds on the event horizon r = rn, we have 

L e ^ = at r = r H . (5.3) 

The solution of the above equation which satisfies the regularity condition on the symmetry 
axis 9 = is 

q? = C(l-cos0), (5.4) 

where C is an integration constant. The above solution implies that the magnetic field 
which can penetrate the event horizon is the only monopole component. By contrast to the 
extreme case, non-monopole components can penetrate the event horizon in the non-extreme 
case, since, in this case, Eq. (15. 2p does not necessarily imply Lq^> = 0. As a result, we can 
conclude that the Meissner-like effect of the extreme black hole appears in the vacuum case. 

In the case that Qp vanishes everywhere, we can obtain global solutions. The solutions 
which satisfy the regularity condition on the symmetry axis 6 = are written in the form 

oo 

= ^2Ri(r)Pl (cos 6) sin 9, (5.5) 

1=0 

where Pi(x) is the associated Legendre function of the first kind with m = 1. Then, the GS 
equation (13. 6p becomes 

a 2 R'l + 2a'aR\ - ^ ^ Ri = 0. (5.6) 

In the extreme case, in order that Ri, R[, and R![ are finite on the event horizon, Ri of / > 1 
vanishes, and the only monopole component / = may remain. 
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B. Of 7^ on the horizon 



In this case, from the regularity condition on the event horizon (14. 3p . we have 

d e ^ = 0, (5.7) 

for both extreme and non-extreme cases. Thus, ^ = is a solution which satisfies the 
regularity condition on the symmetry axis 9 = 0. The magnetic field does not penetrate the 
event horizon at all. 



VI. CASE WITH AN ELECTRIC CURRENT 



In this section, we focus on the extreme case r± = th and assume that the rotational 
velocity of the magnetic field f2p is constant. 

In the case of f2p = 0, it is seen from Eq. (I3.6P that since D = a 2 , the following condition 
should be satisfied on the event horizon: 

dl 2 , . 

i*=°- < 61) 

The above condition allows / = const, on the event horizon. However, since / should vanish 
on the symmetry axis 8 = 0, the allowed constant is zero. Thus, in this case, the same 
argument as was used in the vacuum case discussed in the previous section is also true. The 
allowed configuration of the magnetic field on the horizon is only the monopole component 
(15. 4p . Hence, hereafter, we focus on the case of ^ 0. 



A. Grad-Shafranov equation near the event horizon 

We are interested in the configuration of the magnetic field near the event horizon. In 
order to analyze the GS equation, we introduce the following dimensionless quantities: 



r =: r R y, (6.2) 

* =: r H ^, (6.3) 

e := r H fi F , (6.4) 

=: eJ(^), (6.5) 

2 % =■ e?r?SM. (6.6) 
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Using these quantities, the GS equation (13.61) becomes 

1 T , U W 



where 

V = y[{y-lf-e 2 y^m 2 9], (6.8) 

U = 2(y - 1 - e 2 y 4 sin 2 6)d y i/j, (6.9) 

W = -eV [sin 2^7/; - 5 (^)] . (6.10) 

The regularity conditions on the event horizon (14.51) and (14. 6p become 

Z(V>) + -r- sinedeip = 0, (6.11) 

47T 

(£Z^ 1 

^<VA + — Bwed (d y i(>) = o. (6.i2) 

The regularity conditions on the light surfaces (14. 9p become 

^etp ^— - ey tan 6 sin 6 (1 - ey 2 sinO) d y if) = at y = y LS ±, (6.13) 

sin 26 

where 

Vls± := (6.14) 

Here, we assume that i[) can be written in the form of Taylor series around the event horizon, 
y = 1, as 

oo 

V = ]£V (B) (0)(v-ir- (6-15) 

n.=0 

The coefficients are, in principle, determined by the GS equation (16.71) with the regular- 
ity conditions fIBTTTD . flBTT2l . and (RH3]) . Using the expression (16TT5) . Eqs. fIBTTTD and fl6TT2"|) 
can be rewritten in the forms 

J( ^(o) ) + i_ sin ^^ = o, (6.16) 

^(^ (0) )^ (1) + f sin^ = 0. (6.17) 

If we fix the functional form of we obtain ip(°> and from the above equations, 

and further, we obtain ip^ n > of n > 2 from Eq. (16. 7ft ; in order to get ip( n > for n > 3, we use 



14 



an equation obtained by [n — 2)-times differentiation of Eq. (16. 7p with respect to y. For 
example, for n = 2, by evaluating the GS equation (16. 7p on the event horizon y — 1, we have 

1 dS 



L ^(2) +2 ^(2) + 2cot ^^ (2) 



d9 sin 29 di> 



d9 sin 29 J \sin 2 9 



+ 12e z - e A i) 



(^ (0) )V (2) 



d 2 S 



2 sin 2 # 2 



(0)V,(!) 



(6.18) 



Here, we should again note that X(^>) cannot be freely specified. The functional form of 
must be chosen such that the regularity condition (16 . 1 3[) on the inner light surface is 
satisfied. 



B. Perturbative analysis 

We consider slowly rotating magnetic fields, or in other words, we assume 0<£<1. We 
rewrite the basic equations in the form of the power series with respect to e, and then we 
construct a solution of ip on the horizon, i.e., by perturbative procedures with respect 
to e. Although, as mentioned, it seems to be impossible to determine the functional form of 
I numerically, we can find it by this method. 

In order to construct a perturbative solution for we write 

oo 

^"ty) = j>&V*. (6.19) 

N=0 

Further, we assume 

oo oo 

T(x) = ^l N+l (x)e N and S(x) = ^ S N+2 (x)e N . (6.20) 

N=0 N=0 

From Eq. (14.81) . the location of the inner light surface is written as 

Vls- = 1 + Y J 2n) L (e sin 9) n = 1 + g sin 9 + 2(e sin 9f + ■ ■ ■ . (6.21) 
A — ' n! n + 1)! 

n=l v ' 

Because t/ls- — 1 = 0{ £ )i we can express the quantities on the inner light surface by using 
the quantities on the event horizon. For example, deip at y — y^s- is written as 

B^,0)^ + e^f + ^) + .... (6.22) 
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Since the main purpose of this study is to see the effect of an electric current on the 
configuration of the magnetic field on the event horizon, we focus on ifj(°K For this purpose, 
we rewrite Eqs. ( I6.17P and (I6.18P in more appropriate forms as follows. 

By differentiating Eq. (16.1 6p with respect to 8, we have 



J = 0. (6.23) 



From Eq. (I6.17p . we have 



47r ^(^(o)) + sin ^_^ ln ^(D =0 . (6.24) 
dtp du 



By subtracting Eq. (I6.23P from the above equation, we obtain 



It is easy to integrate the above equation, and we have 

^ = C'sin0^, (6.26) 
do 

where C is an integration constant. In order to obtain ip^\ we use Eq. ( I6.26P rather than 
Eq. (IBTTTD . 

By substituting Eq. ( 16. 16ft into Eq. (16.23}) and using Eqs. (16. 5 p and ( 16. 6p . we have 

The above equation is equivalent to Eq. (I4.4p . By substituting Eq. (I6.27P into the first term 
on the right-hand side of Eq. ( 16. 18ft . we obtain 



L 9 V (0) = e 2 sm 2 



- 12L^ (0) + L g tp {2) + 2tp 



(2) 



n n fdtpW i dS ,, 2 A , m / 1 d 2 S , m 

+2cot# — ^ (2) +ip {1) 2 = -ip {1 > 

V dS sm29dip^ J V V 2sin 2 ## 2 ^ 

We shall use the above equation rather than Eq. ( 16.181) . 



. (6.28) 



1. Zeroth-Order Solutions for ipW 

Here, we obtain the zeroth-order solutions for Hereafter, the arguments of Sn and 
Xtv are tp^ as long as we do not specify them. 
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First of all, we write down the equations to obtain the zeroth-order solutions for ip(°\ 
From the lowest order of Eqs. (16.131) . (I6.16p . and (16.281) . we have 

1 <M (0) 

L e ^ ] = 0, (6.31) 

where from Eqs. (|6"3j) . f lB3|) . and f l6T20|) . 

5 2 = lQn 2 !^. (6.32) 
We can easily integrate Eq. ( 16 . 3 1 f) and obtain 

^ = Cl ) °\l-cos6), (6.33) 
where Cq ^ is an integration constant. The above result implies that ip^ has only the 



1 1 (X) = -^C! ) 0) X{2-X), 



monopole component. Then, substituting Eq. (I6.33j) into Eq. (I6.30p . we have 

1 

Air 
where 

X=^r. (6.35) 



(0)' 

') 

It is non-trivial whether the lowest order of the inner light surface regularity condition 
(I6.29P is satisfied by ip^ and X\ obtained above. From Eqs. (I6.34p and (I6.32p . we have 

S 2 (X) = 16tt 2 J 1 (X)^(X) = 2C$ )) X (2 - 3X + X 2 ) . (6.36) 
dX V / 

It is easy to check that Eqs. (16.331) and (16.361) satisfy Eq. (16.291) . Namely, we have obtained 

a small electric current which satisfies the lowest order of the inner light surface regularity 

condition. 



It is worthwhile to notice the meaning of the zeroth-order solutions for i.e., ip { 



(o) 
• 



In the limit e — > 0, Qp and I become zero, whereas ip^ becomes ipQ . Since the case 
fip = I = corresponds to the vacuum case, i.e., the case without an electric current, ip^ 
corresponds to the vacuum solution. As we showed in §V, the vacuum solution has only 
the monopole component on the event horizon. Eq. (16.331) is consistent with this fact. The 
small electric current / = eX\ can be regarded as a result of the slowly rotating monopole 
field 4 0) - 

17 



2. First-Order Solutions for ip^ 



Next, we consider the correction of 0{e 1 ) to the zeroth-order solution for Hereafter, 
we assume Cj, 0) ? 0. In our perturbative method, the case of Cq ^ = is quite different from 
the case Cf } + 0. If we choose = 0, we can obtain only the trivial solution ip(°> = 
using our perturbative method. We prove this statement in Appendix [Dl 

The equations determining ip^ are derived from Eqs. (I6.13p . (I6.16p . and (I6.28P of 0(e 1 ): 

dip[ 0) 1 ( dS 2 (o) . c 



d9 sin 29 \ dip 

L e ipf ] = 0, (6.39) 



S 3 = 16 7 r 2 (X 1 ^+X 2 ^'). (6.40) 



where 

dip dip 

We can see from Eq. ( 16. 39ft that ipf 1 is also the monopole solution. Thus, the first-order 
correction merely adds a constant of 0{e l ) to the integration constant of the zeroth-order 
solution. As a result, without loss of generality, we may assume for the first-order solutions 
that 

^ 0) = 0. (6.41) 
From Eqs. ( I6.38p . (I6.40p . and the above equation, we have 

X 2 = = S 3 . (6.42) 

We should check the inner light surface regularity condition (I6.13p . Since ip^ = X 2 = S 3 = 0, 
Eq. (I6.37P becomes 

du smzt/ dip 

In order to estimate the above equation, we need ipQ , which is the zeroth-order solution for 
ip m . We obtain ip§ ] from Eq. (KT26]) of O(e ) as 

tfWsin^, (6.44) 
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where we assume that C is the order of unity. Substituting Eq. (I6.33j) into the above 
equation, we obtain 

4 1 ) = C'C^ sin 2 9. (6.45) 

By substituting Eq. (I6.45P into Eq. (16.431) and using the functional form of S 2 given by 
Eq. (I6.36p . we can see that Eq. (I6.43[) is satisfied. Here it is worthwhile to notice that 
ipQ , as well as ip^\ necessarily corresponds to a vacuum solution. It is easy to check that 
Eq. (I6.45P is consistent with Eq. (15. 5p . 



3. Second-Order Solutions for ip^ 

Hereafter, we will make frequent use of Eqs. (I6.29p and (I6.43P without giving an explicit 
reference. Now, we consider the correction of 0(e 2 ) to the zeroth-order solution for ip(°> . The 
equations determining ip 2 can be obtained from Eqs. ( I6.13p . ( 16.161) . and ( 16.281) of 0(e 2 ): 



d9 sin 29 \dip^ 2 



sin 2 fl 



+ 



dip, 



(2) 



2tan# 



1 dS- 



d 2 S- 



2 /,(1) \ ,/,(!) 



2 sin 29 dip' 



2 ,,(2) 



d9 sin 29 dip 



-2tan^o 



(2) 



dli (o) 1 . a dip 2 

V 2 +Is + ^ Sm6 —9 



(o) 



Leipf = sin 2 9 



L 9 4 2) + 24 2) + 2cot^ #0 



(2) 



1 dS 2 ^ 2) 



d9 sin 29 dip 



+ 4 1} 1 



d 2 S. 



2 , (1) 

2sin 2 # dip 2 Wo 



(6.46) 
(6.47) 



(6.48) 



where 



S 4 = lQn 2 { Ii^+X ; 



dX 1 



(6.49) 



dip dip 

In order to derive the above equations, we have used ip± = X 2 = and ipf 1 = obtained 
by substituting ip^ = into Eq. (I6.26p . 

It should be noted that ipQ appears in Eqs. (I6.46P and (I6.48p . In order to determine 
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4 2) , we use Eq. <KT7} of 0{e 2 ): 



L e ^f ] + 2 cot 9 



d4 0) 



i 



d«S 2 (o) c 



d0 sin 20 

By substituting Eq. (I6.46P into the above equation, we have 



0. 



(6.50) 



(0) 



2 cot sin 2 



2tan0 



d 2 S- 



2 sin 20 # 2 



2 ,,(1) \ „/,(!) 



(2) 



1 dS. 



2 ./.(2) 



d0 sin 20 



^-2 tan 0^o 



(2) 



(6.51) 



By subtracting the above equation from Eq. (16.481) . we obtain the equation for ip^ as 



L e ^ ] + 6^ 2) - 24 1] = 0. 



(6.52) 



It is easily seen from Eq. flfT45l) that ip^ = ^?V 2 

is a particular solution for the above equa- 
tion. Thus, the general solution of the above equation is expressed by a linear combination 
of this particular solution and general solutions of the following homogeneous equation: 



L e f + 6f = 0. 



The general solution of Eq. (16. 53ft is given by 

/ = sin0 [c p P 2 1 (cos0) + c q g 2 1 (cos0)] 



(6.53) 



(6.54) 



where c p and c q are arbitrary constants, and P 2 X and Q2 1 are the associated Legendre 
functions of the first and second kinds with / = 2 and m = 1, respectively. From the 
boundary condition at = 0, c q must vanish. Hence, the most general solution of Eq. (16.521) . 
which satisfies the boundary condition at = 0, is given by 



4 2) 



Cf sin 2 cos + ^-^ sin 2 0, 



(6.55) 



where we have used P2 1 = (3/2) sin 20, and is an arbitrary constant. Here it is worth- 
while to notice that the above result also corresponds to a vacuum solution. It is easy to 
check that Eq. (16.550 . as well as if;^ and ipQ , is consistent with Eq. (15. 5p . 
By using Eq. (I6.36P and substituting Eq. (16.450 into Eq. (16.480 . we have 



(0) 



3d 0) C"(l + 2C cos 0) - 4d 2) cos 



sin 4 0. 



(6.56) 



The above equation implies that, in general, the magnetic field on the event horizon includes 
non-monopole components of 0(e 2 ). 
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We can easily integrate this equation and obtain 
sin 2 # 



60 



2 [3C {0) C 2 - cos 9 (3 cos 2 0-7) + 15CfC" (cos 2 9 - 5) 

where we have chosen the integration constant so that ^ |fl=o = ^4 
correction consists of only non-monopole components. 

The functional form of X 3 is determined by using Eq. ( I6.47P as 

1 



(6.57) 
0, i.e., this 



MX) 



x 2 (2-xy 



1207T 

lbC^C ( 



3C^C' 2 



2Q 



(2) 



2 - 18X + 9X' 



(6.58) 



where X is defined by Eq. (I6.35p . 



4- Solution near the inner light surface with corrections up to 0(e 2 ) 



The solution with the corrections up to 0(e 2 ) behaves near the inner light surface as 

i, = 4 0) + e4 0) + e 2 ^ + + e^W) (y - 1) + ^ (y - l) 2 + ■ • • . (6.59) 

Although the solutions for ibff for (n, N) = (0, 0), (0, 1), (0, 2), (1, 0), (1, 1), (2, 0) have been 
derived in the previous sections, we again show them with a slightly different parameteriza- 
tion: 



4 0) 



a 



(0), 
1 



COS I 



0. 



(i) 



,(0) 

C^C sin 2 6, 



1 



Cf j cos 9 (3 cos 2 - 7) + ^C^C (cos 2 0-5) 



sin 2 0, 



(i) 



(2) 



0. 
1 

2 



3 ( C^C' 2 



DCf> ) cos 9 + C^'C 



AO). 



sin 2 0, 



(2) 

where Cg is given in this parameterization as 



a 



(2) 



iod 0) 



By using the same parameterization as the above, the electric current is given by 



X{X) = -—X I 



2C, 



(0) 



(6.60) 
(6.61) 
(6.62) 

(6.63) 
(6.64) 
(6.65) 

(6.66) 



e 2 X (2 - X \ \cfc\\ -X) + 2Cf } (2 - 18X + 9X 2 ) } , (6.67) 
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where X is defined by Eq. ( 16.351) . We see that the arbitrary constants are only C N (N = 0, 2) 
and C . The reason this result takes this form is because if we choose ip and d y ip on the 
event horizon such that the regularity conditions (16 . 1 6[) and (16.171) are satisfied, then ip and 
X are completely determined. 

VII. SUMMARY AND DISCUSSION 

We studied a force-free magnetosphere in a static spherically symmetric black hole space- 
time with a degenerate event horizon. We have found that if an electric current exists, higher 
multipole components of the magnetic field can be superposed upon the monopole compo- 
nent on the event horizon even if the two horizons degenerate into one horizon. This result 
is consistent with the numerical result given by Komissarov and McKinney: they showed 
that the magnetic field lines of higher multipole components can penetrate an extreme Kerr 
black hole if conductivity exists. The detailed geometrical structures of the extreme Kerr 
black hole and the extreme Reissner-Nordstrom black hole are different from each other. 
However, since the degenerate structures of the horizons of these black holes are similar, the 
present results may be applicable to a certain extent for the extreme Kerr black hole. 

If we require that there is no monopole component in the lowest-order configuration on 
the horizon, or equivalently, ^o°^ = 0; we obtain the trivial solution ip^ = 0, even though we 
take all-order corrections into account (see Appendix [D|. Thus, the proposition in Appendix 
D seems to imply that there is no non-trivial configuration without a monopole component 
on the event horizon of the extreme Reissner-Nordstrom black hole, even if an electric current 
exists. But this is not necessarily true. In order to see this fact, note that there is an exact 
monopole solution for the Grad-Shafranov equation ( 13. 6ft : 



where C is an arbitrary constant. By contrast, the proposition in Appendix D implies that, 
if ip^ = 0, there is no higher-order correction by which the configuration of the perturbative 
solution on the event horizon approaches to the monopole configuration in our perturbation 
scheme. In other words, our perturbative solution with vanishing ijj^ cannot approach 



^ = C(l -cos0), 



(7.1) 



with the electric current 




(7.2) 
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to the above exact solution, even if we take into account all-order corrections. This fact 
suggests that even if an exact solution with a non-monopole configuration on the event 
horizon exists, the perturbative solution with vanishing ip^ cannot approach to such a 
solution in our perturbation scheme. This possibility may arise from the assumption for the 
electric current (I6.20p . which may be too strong, though the present analytic perturbation 
studies are impossible without this assumption. 

We would like to stress again that it is very difficult to obtain a stationary force-free 
magnetosphere by solving the Grad-Shafranov equation for the extreme black hole space- 
time numerically. Thus, we need to invoke analytic methods, as in the present study, or 
numerical techniques to follow the dynamical evolution of a force-free Maxwell field until 
a stationary configuration is realized, as Komissarov and McKinney used. As discussed in 
this paper, in the case that there are two light surfaces in the extreme Reissner-Nordstrom 
black hole spacetime, even though the magnetic field is regular, both on the event horizon 
and inner light surface, it will be singular on the outer light surface. If the angular velocity 
of the magnetic field is constant, two light surfaces necessarily exist. Thus, if the dynamical 
evolution of a force-free Maxwell field can be followed until it becomes stationary, then it 
is expected that the angular velocity decays far from the black hole so that the outer light 
surface does not exist. The extremity of charge or angular momentum changes the structure 
of boundary conditions for the Grad-Shafranov equation and seems to strongly affect global 
structures of the black hole magnetosphere. 

Finally, we would like to suggest that analytic solutions obtained by this perturbation 
scheme becomes a benchmark for a numerical scheme to obtain solutions for stationary con- 
figurations of astrophysical magnetospheres, since our perturbation scheme is also suitable 
for non-extreme black hole cases. 
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Appendix A: Derivation of Grad-Shafranov equation 



A stationary, axisymmetric electromagnetic field implies dtA a = = d^Aa. Then, from 
Eq. (13.31) and the force- free condition (13.4)) . we have F te /F eip = F tr /F rif . Using these equa- 
tions, the components of F a b in the static coordinate system (12 .ip are written in the form 

/ tt F d r A v tt F d e A v \ 

-d r Ap -d e A v 

-Q F d r A v d r A v ^B* 
\ -n F doAy d g A v -^jB* 

where fl F is defined by Eq. ( 13. 5p . and 

r 6 sin 2 9 



(Al) 



/ 



7 := 



(A2) 



is the determinant of the intrinsic metric of the spacelike hypersurface labeled by t. 

Note that Q F can be regarded as the angular velocity of the magnetic field. We consider 
an observer with an angular velocity dtp/dt = Q F . His or her 4- velocity is given by u M = 
T(l, Q F , 0, 0), where T is a normalization factor. The electric field for this observer is given 
by E a = F a bU a , and we can easily see from Eq. ( jAlf) that E a vanishes. Thus we may say that 
this observer is co-moving with the magnetic field, and the angular velocity of the magnetic 
field is VL F . 

Substituting Eq. (IAip into the Jacobi identity d[ a F hc ^ = 0, we have 



(d r n F )d e A v - (d 9 n v )drAp = o. 



(A3) 



The above equation implies that d a Q F oc d a A^, or equivalently, the equi-^F surface agrees 
with the equi-A^ surface. Thus we have 



n F = Q F {A V ). 



(A4) 



Using Eq. (IAip . the Maxwell equations imply the following equations: the t-component 
implies 

d r (r 2 Q F sin 6d r A v ) + d e r d e A v J = -Ana^J*; (A5) 

the yj-component implies 



( -?^d r A v ) + d e ( 1 ^— Q deA v ) = -Ana^J^ 
sin v I \ A sin V 



(A6) 
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the r-component implies 

d e (aB v ) = 4 7 ra v ^J r ; (A7) 

and the ^-component implies 

dr(aB v ) = -4vra v ^J 9 , (A8) 

where 

B ip = r 2 sin 2 6B*. (A9) 
Using Eq. ( lAljl . the force-free condition implies 

J r d r A v + J e d e A v = 0, (A10) 
{J* - JtQ^drAp + ^jB*J e = 0, (All) 
( J v - j'Q^deA^ - y/yBT = 0. (A12) 

Substituting Eq. (JSHJl to Eq. flATll . and substituting Eq. (JSID to Eq. (TAT21) . we have 

(J* - J*ft F )<9 r A, - i — B^driaBp) = 0, (A13) 

(J v - J'fip)^ - -^-^^(aS^) = 0. (A14) 
Ana 

From the above equations, we have 

d r {aB v )d e A v - d e (aB v )d r A v = 0. (A15) 

The above equations imply 

aB v = B{A V ). (A16) 
Using the above equation, Eqs. ( 1A13I) and ( 1A14I) imply 



1 ^7 

- J'fip = -4 T-rP-^- ( A1 7) 

47ra 2 r 2 sin 2 # cL4^ 



From Eqs. (IA5j) and (1A6I) . we have 



Or . n + -rd e . . y H t A(a r A)9 r fi F + (<%A)<9 e fi F 

\ smv J A \ sma / A 

= -Ana^ {J v - J l n F ) , (A18) 

where D is defined by Eq. (13. 101) . Noting that Qp is a function of A^ and substituting 
Eq. (1A17I) into the right hand side of Eq. ( 1A18j) . we have 



d 2 .A v + 1 + ^ ) = 0, (A19) 
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Polar Cap 



r(l)cosd(l) 
ft cos ft 




r(t)s'md(l) rcsinft 



FIG. 3: The schematic diagram of a polar cap which is an axisymmetric two-dimensional spacelike 
surface parameterized by the proper length t measured from the symmetry axis 8 = along the 
polar cap. 



where 



and 



LtiA 



sin 6d e ^—jr 



r 2 dVl 2 

N := A(d r A v )d r D + (dgA^dgD + - sin 2 6-^ [A(0 r A„) 2 + (d e A^) 2 } 



r 2 dB 2 
~2~dA~ n 



(A20) 



• (A21) 



The above equation is called the Grad-Shafranov equation. 



Appendix B: Electric current and magnetic flux 

Here we introduce electric current I and magnetic flux $ on a spacelike hypersurface 
labeled by t which penetrate downward and upward an axisymmetricpolar cap, respectively. 



These quantities were first introduced by Macdonald and Thorne js| and are related to B 
and as follows. The polar cap is parameterized by I and <p, where I is the proper length 
on the polar cap from 6 = 0. The coordinates r and 9 on the polar cap are given as functions 
of £, i.e., r = r(£) and 9 = 9(£): by definition, 9(0) = 0, and we assume that r(0) > r + (see 
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Fig. [3]). The orthonormal tangent vectors of the polar cap are 

dr dQ s 



dt dt ) ' 
0, 



r sm i 



(Bl) 
(B2) 



Then, the upward unit normal to the polar cap is 

r 2 ( _ dO dr\ 



0. 



- 



di 



(B3) 



We assume that the edge of the polar cap is r = r e and 9 = 8 e . Then, denoting the 
proper length £ at the edge by £ e , we have 

»2tt 



JO 



« J*n,r sin 



dd ^ . . dr 
—d e (aB v ) + — d r {aB v ) 



dt 



d(aB^ 
di 



-aBu 



(B4) 



(r,8)=(r c ,e c ) 

where we have used Eqs. flA7j) and (IA80 in the second equality and assumed that B v \ e=0 = 
from the regularity requirement. We can easily see from the above equation that the electric 
current through the polar cap is a function of the coordinate values of the edge, (r e , 8 e ). By 
a similar consideration to that for the electric current /, the magnetic flux \1/ can be written 
in the form 



2tt 
</0 



s ljk F jk riir sin 9d£d(f 



2tt 



2tt 



d6 



dr 



2tt 



o 

lc dA v 



dd 



-F, 



ipr 



dr \ r 3 sin 9 



di 



d£d(p 



di = 2-kA 



o 



<p\(r,8)={r e< e e ), 



(B5) 



where e^ k [e we = 1/ y/j) is the components of the skew tensor in the spacelike hypersurface 
labeled by t, and we have assumed A^l^o = from a regularity requirement. Thus we have 



/ = — B and 
2 



* = 27rA v . 



(B6) 



Rewriting Eq. ( 1A19I) using I and the Grad-Shafranov equation Eq. ( 13.61) is obtained. 



Appendix C: Relation between the Kerr-Schild and Static coordinate systems 



The line element of the Reisner-Nordstrom spacetime with the Kerr-Schild coordinate 
system (T, <3>, R, 0) is given by 
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R 2 



R 2 + 2MR - Q 



-dT 2 + R sin 2 Qd& 



2MR - Q 2 



-AT 



R 2 d& 2 



(CI) 



R 2 V R 2 + 2MR-Q 2 

where M = r + + r_ and Q 2 = r + r^. The relation between the Kerr-Schild coordinate system 
and the static one is given by 



dT 

<i$ = dip, 

dR = dr, 

dQ = dO. 



dt -\ : — — c/r, 



(C2) 
(C3) 
(C4) 
(C5) 



From the above relation, we have $ = ip, R = r, and 6 = 9, and 



d_ 

&T 
d_ 

a¥ 

d_ 

dR 
d_ 

dQ 

Using the above relations, we have 



d_ 

di' 
d_ 

dip' 

2Mr - Q 2 d 



A 



d_ 
d9' 



A^ — A§. 







dt dr ' 



(C6) 
(C7) 

(C8) 
(C9) 

(CIO) 



By virtue of the stationary, axisymmetric nature of the electromagnetic field, we can easily 
see that the components of F a b in the Kerr-Schild coordinate system are given as 



Ft<s> 
Feb 



0, F TR = £l ¥ d r A v , F Te = n F d e A^, 

F<s>e = —df)A 



-d r A^p, 



- -(2Mr - Q 2 )n F d e A^ 
1 



Asin6> 



B - —(2Mr - Q 2 )tt F sm9d e A v 



(Cll) 
(C12) 

(C13) 
(C14) 



It should be noted that all the ordinary derivatives of the Kerr-Schild coordinates are 
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equivalent to those of the static coordinates for the stationary axisymmetric field A$ = A^, 
d T A<s, = dtAp = 0, dq,Aq, = d v A v = 0, 

d R A* = d r A v , d A^ = d A v , d 2 R A$ = d 2 A 9 , and <9|A $ = d 2 A v . (C15) 

Thus, even if the Kerr-Schild coordinate system is adopted, the equation for A® takes exactly 
the same form as Eq. (1A19|) . 



Appendix D: Non-existence of non-monopole solution 

Proposition: Within the perturbation scheme developed in this paper, the solution for 
with vanishing lowest-order solution i/j^ is the only trivial solution ip^ = 0. 

Proof. We prove this by induction. We have already shown that, if we impose the non- 
existence of a monopole component for the lowest order of the perturbative solution, we 
obtain 

4 0) = 0. (Dl) 
Here, we assume that tpff = for < N < M, or equivalently, 

oo 

V^0)=s m+1 X>^at + m + i(0)- (D2) 

N=0 

Then, we have 

oo / oo \ 

Jty<°>) = £ e N X N e M ^ Yl £ ^m + i = e M+ %'^S +1 + 0{e M+ \ (D3) 

N=0 V J=0 / 

where 



x :- 



(D4) 

and we have used Xtv(0) = 0, which is required from the regularity condition at 9 = 0. 
Substituting the above equation into Eq. (16.161) . we obtain an equation of 0(e M+1 ), 



t'/( ) 1 ■ n ?t\k.\ 

X ^V+i + ^ sm 9 dQ = 0. (D5) 
Integrating the above equation, we have 

r S u,2fl i~ 2 ^' 

(D6) 



?/ .(o) _ r (o) 



'1 + COS0) 2 
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where C^ +1 is an integration constant. The regularity condition implies ip$ +1 \e=o = = 
dij>M+i/ d&\e=o- Hence, if Cj^+i does not vanish, 

- 2nl ' = 1 + c 2 (D7) 

must be satisfied, where c is an arbitrary constant. 

In the neighborhood of = 7r, sin ~ n — 8 and cos ~ -1 + (0 - tt) 2 /2. Hence, if c£ } + i 
does not vanish, we have 

lim^+^oo. (D8) 

If we require the finiteness of ^m+i at = 7r, then Cj^+i must vanish, and, as a result, 
■0M+1 = is obtained. Q.E.D. 
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